function b=SumPow(a,n,dim)
if nargin==1
    n=2;
    dim=2;
elseif nargin==2
    dim=2;
end
b=sum(power(a,n),dim);
